# usage: python /home/orsi/brahmsProject/analysisCode/calcSpontCurvFromLppMomentsLog.py

import sys,string,math

# Avg and SE for vals from force.dat
tau1m = []; # lists to collect first integral moments of monolayer
f = open('lppMoments.log','r')
while 1: # Create infinite loop
        line = f.readline()
	if line.find(' TORQUE TENSION ') != -1: # line contains P1
		words = string.split(line) # Split each line into columns
		print words[7]
		tau1m.append(float(words[7]))

print len(tau1m)
